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p ^ 1 [ Abstract 

\ We introduce a notion of viscosity solutions for the two-phase Stefan problem, which 

■ incorporates possible existence of a mushy region generated by the initial data. We 

show that a comparison principle holds between viscosity solutions, and investigate the 
coincidence of the viscosity solutions and the weak solutions defined via integration by 
parts. In particular, in the absence of initial mushy region, viscosity solution is the 



unique weak solution with the same boundary data. 







> 

in 

^ ; 1 Introduction 

The classical two-phase Stefan problem models the evolution of temperature in a material 
with two distinct phases (solid and liquid). The solid-liquid phase transition occurs at a 
given constant temperature, which we set at zero. The problem is described by an enthalpy 
function h, which represents the internal energy density in each phases (see Meirmanov 
[19] and Oleinik et al. [20]). More precisely, let H C M" be a bounded domain in R n 
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with a C 2 boundary. Then h(x,t) : Q x [0,T) 



solves the following problem: 



fH = A X (h) 



in n x (0,T), 



h(x, 0) = /io(s) in ft, 



(ST) 



Ufa) 



on dtt x (0,T). 



In this paper we will set 6(x, t) as a continuous function on d£l x [0, T] and ho as a bounded 
measurable function with x(^o) continuous on 0, such that x(7io) = G(-,0) on 50. Here 
u(x,t) := x{h){x,t) describes the evolving temperature, where x is defined by 



X(h) 



h. 

0, 

h+1, 



h>0, 

-l<h<0, 
h < -1. 



(1.1) 



Note that ft, is determined by it almost everywhere if T(u) := {u = 0} = {—1 < h < 0} 
is of measure zero, but not in general. The interior of zero temperature set T(u), if exits, 
is called a mushy region. It is well known that mushy region might appear in a Stefan 
problem with heat sources, see Bertsch et al. [4J. When there is no heat source, the 
mushy region can only be generated from the initial data. Indeed if |{ — 1 < ho < 0}| = 
(here | • | denotes the Lebesgue measure in IR n ), it is shown in [21] (see also [16] and [1] ) that 
|{ — 1 < h(-,t) < 0}| = for all t > and ()STj) can be written in terms of the temperature: 



u t 



Au = 



in {u > 0} U {u < 0}, 



V n = \Du + \ - \Du~\ on r(it) = d{u > 0} = d{u < 0}, 

u(x,0) = u (x) := x(ho(x)) 

u = 9 on dnx(0,T) 



(STn) 



Here u + and u denotes respectively the limit taken from {u > 0} and {u < 0}, and V n is 
the outward normal velocity of the free boundary with respect to {u > 0}. 

On the other hand, if |{— 1 < ho < 0}| > with no heat source then the mushy region is 
non- increasing as the solution evolves (see [21]): in this case the zero temperature set T(u) 
may exhibit instability such as infinite speed of propagation as its interior diminishes (see 
[5] where the interface evolution is investigated in the case of one-phase Hele-Shaw flow 
with initial mushy region). 



For both ()STp and (jST-up . even solutions with initially smooth configuration may develop 
singularities in finite time, due to merging and splitting of the liquid-solid regions. Thus it 
is necessary to introduce a generalized notion of global-in-time solutions. 
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The first global notion of solutions for (|ST|) has been derived by a weak formulation via 
integration by parts: see Ladyzenskaja et. al. [18] or Friedman [15]. This approach 
allows the presence of initial mushy region. It was shown there that a unique weak solution 
exists for bounded measurable ho and 9. Moreover the temperature u(x,t) turns out to 
be continuous when ho and 9 are sufficiently smooth - see Caffarelli & Evans [H] and 
DiBenedetto [13] . 

More recently, the notion of viscosity solutions, first introduced by Crandall & Lions 
[12] , has been proved useful in the qualitative study of free boundary problems, due to 
its flexibility with nonlinear PDEs and free boundary conditions: see e.g., [2]-|3] and [7] 
which address a general class of free boundary problems including (jSTttj) . and also see |14] 
where the heat operator is replaced by a class of uniformly parabolic operators. However, in 
the two-phase case, the equivalence between viscosity solutions and weak solutions, and in 
particular the uniqueness of the viscosity solutions of (jSTj) . has not yet been investigated, 
even in the case of (jSTnj) . (The only notable instance is [9] where the authors show a 
comparison principle between viscosity subsolutions and more regular supersolutions, so- 
called "R-supersolutions" . Even thought their result could be used to simplify the proof of 
Lemma 13.81 it is restricted to the heat operator, whereas our method is more robust and 
extends to a wider class of problems, see Remark 11.21 ) In this paper we address precisely 
this point. We introduce a notion of viscosity solutions for (|ST|) . to derive the following 
theorem: 

Theorem 1.1. (a) There exists a viscosity solution of (jSTj) in Q. Furthermore, a com- 
parison principle holds between viscosity solutions and thus maximal and minimal 
viscosity solutions exist. They respectively correspond to weak solutions of (ISTj) with 
maximal and minimal initial enthalpy for given temperature uq . 

(b) Any weak solution of (ISTj) in the sense of [8] is also a viscosity solution. 

(c) If |{ — 1 < /io < 0}| = 0, then there is a unique viscosity solution of (ISTuj) which 
coincides with the weak solution of (jSTj) . and also coincides with the notion of viscosity 
solutions discussed in [12] and [2j. 

Previously, one of the authors investigated in [17] viscosity solutions theory for the one- 
phase version of (jST^P (also see [10] and [6] for the case of degenerate diffusion). While 
one-phase Stefan problem is a special case of the two-phase Stefan problem, this work does 
not immediately apply to the two-phase version: we will point out the differences as we 
proceed in the proof. For example, the free boundary evolution for the two-phase setting is 
no longer monotone (that is, the positive set of u may shrink or expand over time), which 
requires a more careful construction of test functions. Moreover, since we are allowing the 
initial data to have mushy region, T(n) may have positive measure with infinite speed of 
propagation, and thus must be treated differently. Indeed, to the best of our knowledge, 
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Theorem 11.11 is among the first results which addresses uniqueness of viscosity solutions for 
a two-phase free boundary problem. 



Below we give an outline of the paper: 

In sections [2] and [3] we introduce a definition of viscosity solutions for (|STp , and prove that 
viscosity solutions satisfy a comparison principle for strictly separated initial data. Note 
that this does not immediately yield uniqueness except under some geometric conditions on 
the initial data (see Corollary I3.9p . however it generates a maximal and minimal viscosity 
solutions (see Proposition 14. 4p . These maximal and minimal solutions indeed correspond 
to the weak solutions with maximal and minimal initial enthalpy. 

In section HI we show that weak solutions are viscosity solutions. Using this result, existence 
of viscosity solutions is established. Furthermore, using the uniqueness of the weak solutions 
and barrier arguments near the parabolic boundary of Q (see Lemma l4.5p . we establish the 
uniqueness of viscosity solutions when \{uq = 0}| =0 and 9 < 0. 

Remark 1.2. The proof of the comparison principle (Theorem \3.1\) relies only on the 
comparison principle available for the parabolic operator ut — Au and monotonicity of the 
boundary condition. It therefore applies to the generalized version of (jSTup described in 
[2], where the second equation in (jSTnP is replaced by a general free boundary velocity of 
the type 

V = G(x,\Du + \,\Du-\,ri), 
where r/ is the spatial (outward) unit normal vector ofTiu), G(x,a,b,r]) is continuous and 

— -— > c 
da ' db ~ 

And in fact, a general uniformly parabolic operator 

u t — F(t, x, u, Du, D 2 u) 

discussed in [11] can be considered in {u ^ 0}. That is, we require that F is an uniformly 
elliptic operator such that 

X\\N\\ < F(t, x, v,p, M + N) — F(t, x, v,p, M) < K\\N\\ 

where t e [0,T], x e O, v G R, p e W 1 , M,N G S(n), N > 0, and < A < A. The 
use of Hopf's lemma at the end of the proof of Theorem \3.1\ can be omitted if we consider 
the inf- convolution W in (|3.ip on shrinking balls as in jl7| . which only requires a more 
demanding notation. 

As for existence, a slight modification of Perron's method presented in |17] combined with 
Theorem \3.1\ yields the existence of a viscosity solution for the generalized version. This 
approach is not taken here since we are also interested in the coincidence of weak and 
viscosity solutions. 
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2 The Stefan problem 



Let &lQ = dfl x [0, T] be the lateral boundary of Q and dpQ = (Q x {0}) U <9lQ be the 
parabolic boundary of Q. 

We shall use two notions of generalized solutions for (ISTj) : weak solutions and viscosity 
solutions. 

Definition 2.1. A bounded measurable function h is called the weak solution of (|STp with 
initial data ho £ L°°(f2) and boundary data 6 £ L°°(dLQ) if it satisfies (see |19j ) 

/ (% + x(/i)A^)) dxcft- / 6^f-dsdt+ [ h (x)ip(x,0) dx = (2.1) 
/or aZZ £ W^' (Q) vanishing on 8lQ and at t = T. 

Due to [15j[T8l[19] ) there exists a unique weak solution /i of problem (|ST|) . defined in Defi- 
nition [2J3 with initial data ho and boundary data 6. Furthermore a comparison principle 
holds between weak solutions ([H]): see the proof of Lemma [4. II for more precise statement. 

As mentioned in the introduction, the temperature u(x,t) is related to the enthalpy h 
through u = x(h)- Formal computations yield that (jSTp provides the free boundary velocity 

V n = \Du + \ - \DvT\ on d{u > 0}nd{u < 0}; 

V n > \Du+\ on d{u > 0} \ d{u < 0}; (V) 

V n <-\Du~\ on d{u < 0} \d{u > 0}, 

where V n is the outward normal velocity of the free boundary of {u > 0}. The inequalities 
are due to the (possible) presence of the mushy region and the non-constant distribution of 
h(x,t). Definitions 12.41 and 12.51 incorporate (jV]). 

Definition 2.2. Set D C BL n+1 is a parabolic neighborhood of a point (xo>*o) */ there is 
an open set U , {xo,to) € U, such that D = U H {t < to} and D € C 0,1 , i.e. the boundary of 
D is locally a graph of a Lipschitz function in x and t. 

Definition 2.3. A continuous function <j) is a classical subsolution (in a parabolic neigh- 
borhood D) of flST]) if 

(a) <f> is Clj in DD {4> > 0} and in D n {<p < 0}, 

(b) r(0) = a{0 > o} = a{0 < o}, 
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(c) \D<f> + \ and \D<f> \ are nonzero on DnT((j)) (this makes T*((j>) smooth), 

(d) <f> t - A(f> < in D n {(/) / 0}, 

(e) the outward normal velocity V n of D n T((f)) with respect to {<j> > 0} is less than or 
equal to \D(p + \ — \D<p~\. 

Similarly one can define a classical supersolution in a parabolic neighborhood. 

Definition 2.4. (a) A lower semi- continuous function v is a viscosity supersolution of 
(ISTp in Q if (i) v(x, 0) > x{ho), (ii) v > 9 on diQ, and if (Hi) for any classical sub- 
solution 4>(x,t) in a parabolic neighborhood D C Q such that <fi < v on the parabolic 
boundary of D, (f><v in D. 

(b) An upper semi- continuous function u is a viscosity subsolution of (jSTp in Q if (i) 
u(x,0) < x(ho), (H) u < 6 on 8lQ, and if (Hi) for any classical supersolution <fi(x,t) 
in a parabolic neighborhood D C Q such that <f> > u on the parabolic boundary of D, 
<J) > u in D. 

Observe that if u is a viscosity subsolution of (]STp . then — u is a supersolution of (jSTp . 

Definition 2.5. A continuous function u is a viscosity solution of (|ST|) if it is both viscosity 
subsolution and supersolution. 

Remark 2.6. Note that above definition of viscosity solutions, unlike that of weak solu- 
tions, does not take into account the enthalpy function h(x,t). In particular, we cannot 
expect uniqueness when the initial data uq satisfies \{u$ = 0}| > 0. A modified definition, 
incorporating the presence of h(x,t) and therefore yielding a unique representation of the 
interface evolution, is under investigation by the authors. 

3 Comparison for strictly separated initial data 

A central property for the viscosity solution theory is the comparison principle, which 
is stated in the theorem below. We mention that a corresponding result holds for weak 
solutions of (jSTp (see [IS])- As mentioned in the introduction (see Remark 1 1.2p . the proof 
presented in this section extends to a general class two-phase free boundary problems with 
nonlinear free boundary velocity. 

Theorem 3.1. Let u and v be respectively a viscosity subsolution and supersolution on Q 
such that u < v on dpQ. Then u < v in Q. 
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First we introduce several notations. Let B r (x,t) denote the closed space-time ball with 
radius r and center at (x,t), 

B r (x,t) = {(y,s) : \(y, s) - (x,t)\ < r}, 

and D r (x,t) the closed space disk with radius r and center at (x,t), 

D r (x, t) = {(y, t) :\y-x\< r}. 

We also speak of the sets 

dD r (x,t) = {(y,t) : \y - x\ = r} and D°(x,t) = {(y,t) : \y - x\ < r} 
as the boundary and the interior of D r (x,t), respectively. 
For any r > 0, define the sup- and inf-convolutions (see |17j ) 
Z(x,t) = sup U(y,s), U(x,t) 

B r {x,t) 

W{x,t)= inf V(y,s), V(x,t) 

B r (x,t) 

For any p > 0, define 

Q p = {x £ Q : dist(x, dfl) > p}, 
The convolutions (|3. If) are well defined in C r , 

C r = Q 2r x (r,T - r). 

Moreover, we fix r > small enough so that Z < W on dpC r . Note that Z, U, resp. W, V 
are also viscosity subsolutions, resp. viscosity supersolutions of (|ST[) in C r . 

It will be convenient to define $ r (x,t) to be the closed set 

* P (M)= (J ^(y,*)- (3-2) 

(j,s)eB r (i,i) 

Lastly, we need the following definition: 

Definition 3.2. Let E be a closed set, let P be a point on the boundary of E, P £ and 
/ei B be a closed ball containing P. The ball B will be called an interior ball of E at P if 
B C E, and it will be called an exterior ball of E at P if B° Pi E = 0. 

Observe that, by definition, the sets {Z > 0} and {W < 0} have interior balls of radius r 
at each point of their boundaries. These balls are centered on the boundaries of {U > 0} 



sup u(y,s), 

D r (x,t) 

inf v(y,s). 

D r (x,t) 



(3.1) 
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and {V < 0}, respectively. This property of Z and W will be essential in the proofs in this 
section. 

To prove that u and v stay ordered, we will argue by contradiction: if u > v at some time 
to, it must be true for Z and W even earlier. Suppose that there is a time to when the 
functions Z and W first hit, 

t = sup {t > r : Z(>, s) < W(-, s) for all r < s < t}. 



Denote Pq = (xo,to) the point of maximum of Z — W at t = to- Note that Z(xo,to) > 
W(xQ,to) because Z — W is upper semi-continuous. 

Note that the positive set of Z (and the negative set of W) may expand discontinuously in 
time, in the event that the zero set of Z (and W) is of positive measure: this is because 
our definition only considers test functions which are nondegenerate in both positive and 
negative phases. Therefore one should put extra care in describing the location of Pq and 
the geometry of the free boundaries nearby. 

Lemma 3.3. The set {u > 0} cannot expand discontinuously, i.e. 

{u > 0, t = T} C {u > 0, t < T} VT > 0. 
Analogous results hold for {v < 0}, {Z > 0}, {U > 0}, {W < 0} and {V < 0}. 



Proof. Suppose that for some T > there is a point P € {u > 0, t = T} such that P ^ 
{u > 0, t < T}. In other words, there exists 5 > such that JB 5 (P) n {ti > 0, t < T} = 0. 
For simplicity we will translate the coordinates so that P = (0,0). Define M = hi&xqu. 
Note that the ball B$(P) contains a cylinder 

C={(x,t):\ X \<h + i h -|-n(^,£)<t<0} 

with h and r\ small positive constants. The situation is shown in Figure [TJ We will construct 
a classical supersolution of (1ST|) in C that crosses u from above to obtain a contradiction. 

Define uj = 2 max j^J- Consider a radially symmetric function in C defined by the 
formula 

'a(|x| 2 — (h — u)t) 2 ) for \x\ < h — cut, 



<j)(x,t) 

' 2M 

V 

where a > is a small constant specified below. 



( (\x\ — h + ut) for \x\ > h — ut, 
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Figure 1: Situation at point P 



Then </> satisfies 

' 2M ( , , 

M 



' 2M ( u _ n-1 j > q { OT \x\>h- (tit, 



(dt-A), 



2a ((h — ujt)uj — n) > for \x\ < h — tot. 



Next, note that the free boundary T(<fi) at time t is the sphere with radius h — ut, shrinking 
with normal velocity oj. Therefore we have 

V n = uj > — - \D6~\ = \D(f> + \ - \D(f>~\ . 

V 

Hence <fi is a classical supersolution of (1ST|) for any a > 0. 

Lastly, note that the initial boundary D = C n {t = —rj/u}} of the cylinder C is a closed 
subset of {u < 0} and therefore maxpu < 0. We can hence take a > small enough to 
make <f> > u on D. 

Thus we have constructed a classical supersolution on C such that <^> > u on the parabolic 
boundary of C, while </> < < u at P = (0, 0). That shows that u crosses 4> from below 
inside C, which contradicts the definition of u. 

The conclusion for {v < 0} = {— v > 0} can be obtained by repeating the proof for —v. 
Corresponding arguments apply to the sets {Z > 0}, {U > 0}, {W < 0} and {V < 0} as 
well. □ 

Lemma 3.4. The point Pq lies on the intersection of boundaries of {Z > 0} and {W < 0}, 

Po e d{z> o}nd{w < o}, 

and 

W(P ) <o< z(p ). 



t 




^Z>0 



\ 



\V <0 



B 2 



x 



Bo 



Figure 2 



Situation at Pq 



Proof. We know that W(Pq) < Z(Pq) from the definition of Po. There are three cases: 

(a) < W(-Po) < -^(-fb): in this case > in some neighborhood of Po ({W > 0} is 
open as W is a lower-semicontinuous function) and W is a supercaloric function in 
this neighborhood; this contradicts the fact that Z cannot be crossed from above by 
any positive caloric function (f> (see 

(b) W(Pq) < Z(Pq) < : the same as above for Z subcaloric in a neighborhood of Po. 
Therefore we arrive at 



Inequality (|3.3p immediately yields that Pq G {Z > 0} n {W < 0}. Now suppose that Po 
lies in the interior of the set {Z > 0}. That implies Pq ^ {W <0,t< t^} since Z < W 
for t < to- But that is a contradiction with continuous expansion of {W < 0}, Lemma 13.31 
The same happens when Pq lies in the interior of {W < 0}. Hence the conclusion of the 



The lemma states that the situation at point Po looks like the one shown in Figure [2j More 
precisely, there exist P\ = (xi,t\) £ d{U > 0} and P^ = (^2^2) G d{V < 0} such that 
[/(Pi) = Z(P ) and V(P 2 ) = W(P ). At P , Pi = P r (xi,ti) = P r (Pi) is an interior ball of 
{Z > 0} and B 2 = B r (x 2 ,t 2 ) = B T (P 2 ) is an interior ball of {W < 0}. Moreover, the ball 
B = B r (x ,t ) = B r (P ) is an exterior ball of {U > 0} at Pi and {V < 0} at P 2 . 

Since Z < W for t < to, the space projections of Pi and P2 are tangent: 



W(P ) < < Z(P ). 



(3.3) 



lemma. 



□ 



B 1 r\B 2 n{t = t } = {P }. 
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T(4>) 



Figure 3: Local situation at P[ 



Therefore the space projections of Pq, Pi and P2 are on a line and we can change the 
coordinates so that 



Lemma 3.5. Suppose that Z(Pq) = 0. Then (xi,t\) ^ (xo,to + r). In other words, the line 
(Pq,Pi) cannot be vertical (i.e., be parallel to the time- direction) and in the future of Pq. 

Proof. Suppose that (P\,Po) is vertical, i.e. x\ = xo, and t\ = to + r - Thus there is 
P[ = (x'y,ti) £ D r (Pi) such that u{P[) = 0. Lemma 13.41 implies that Pq lies on the 
boundary of {Z > 0}. Therefore u < in the interior of $ r (Po) ? defined in (|3.2j) . This and 
Lemma [3.3l vield that u < in the interior of D r {P\) and P[ lies on the boundary of D r {P\), 
i.e. \x\ — x']\ = r. We will construct a classical supersolution in a parabolic neighborhood 
of P[ that crosses u to get a contradiction. 

First, let uj = \fr/3 and h = 1 + co and consider the set 

E := {(x,t) : -1 < t < 0, 1 - w - < |x| < 1 + 2uj - ujt}. 

(see Figure [3]). 

We only show the construction for n > 3. The computation is similar when n = 2. Thus 
consider a radially symmetric function (j)(x,t) on E of the form 



P = (0,0), Pi = (diei.ti -t ) and P 2 = (-d 2 e u t 2 - t ). 



(3.4) 




a{— \x 
b(—\x\ 



2-n 



+ (h — ojt) 2 n ) when \x\ > h — cot 
+ (h — ojt) 2 ~ n ) when \x\ < h — ujt 



where a = 



2n-4 



h n and b > is a small constant to be chosen later. 



then satisfies 



(0 t - A) <t> = 8 t (j) = c(n - 2)u(h - ujt) 1 - 11 > 



if \x\ ^ h — ojt 
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where c = a or b. The free boundary of eft is the shrinking sphere r(</>) = {\x\ = h — cot} 
and therefore the free boundary condition reads as 

V n = u> a(n - 2){h - ut) 1 -" 1 - \D<j)~\ = \D(j) + \ - \D<j)~\. 

Hence is a classical supersolution on E. Now define 

5 = min (j> > 0. 

\x\=l+2u-cut 
-1<*<0 

and choose e, < e < r, such that u < <5 on B4 £ (P[). This can be done since u(P[) = 
and u is upper semi-continuous. Let E £ and (f) e be rescaled £7 and (f) under the parabolic 
scaling (x,t) = (ex + ex\jr + (1 - e/r)x' 1 ,e 2 i + t\). Note that £7 £ C B^(P{) and (/> e is still 
a supersolution in E 12 . 

Finally, observe that the distance of the lateral boundary of <& r (Po) from the time axis 
{x = 0} in the rescaled coordinates can be expressed as 

|x| = V-2ri-e 2 i 2 + 1 > \f^ri+ 1 if - 1 < i < 0. 

Thus clearly u<0on||x|<l + 2w, t = — 1| and so we can choose 6 > small enough so 
that 4> £ (x,t) > u(x,t) on the part of parabolic boundary of E e where 4> e < 0. This can be 
done as this is a closed subset of {u < 0} and u is upper semi-continuous. 

With above choices of b and e, (fr £ is a classical supersolution of (|ST|) such that <f> £ > u on 
the parabolic boundary of i? 5 and <j) e (P[) < = u(P[). Therefore u — cj) e has a maximum 
inside E e , which yields a contradiction with u being a viscosity subsolution. 

The proof can be repeated with — W in place of Z and other appropriate substitutions. □ 

With the previous results, it is now possible to establish the following lemma: 
Lemma 3.6. 

Z(P ) = W(P ) = 0. 

Proof. It is enough to prove Z(Pq) = because W(P$) = then follows by repeating the 
proof for —W. 

Due to Lemma 13.41 Z(Pq) > 0. Suppose Z(P$) = L > 0. As in the proof of Lemma 13.51 
there is a point P[ G D r (P\) such that u(P[) = L. First we discuss a couple of extreme 
cases: 

(a) If (Pi,Pq) is vertical with Pi in the past of Pq, arguing as in the proof of Corollary 
13.71 we see that (P2,Pq) must be vertical with P2 in the future. So either W(Pq) = 0, 
but that cannot happen by Lemma 13.51 or W(Pq) < 0. For W(Pq) < 0, one could 
argue as in (b) with — W in the place of Z (and other appropriate substitutions). 
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Figure 4: Construction of a supersolution when Z(Pq) > 0. 



(b) If (Pi, Po) is vertical with Pi in the future, the fact that P$ S d{Z > 0} implies u < 
in the interior of <£ r (Po), therefore u < in the interior of D r {P\) by Lemma 13.31 and 
P[ must be on the boundary of D r (P\). The rest of the argument is parallel to the 
one shown below. 

Now suppose that (p},Pi) is not vertical and t\ < tQ. This allows us to find a small 
r > and points Q = (x,ti) and Q T = (x,ii — r) such that Df l (Q) n {u > 0} = {Pi} and 
Dh T (Qr)r]{u > 0} = 0, where h = \P\ - Q\, < h - h T « 1 (see FigureHD. For simplicity 
we translate coordinates so that Q = (0, 0). 

Denote M = maxg u. We will now construct a positive classical supersolution (ft of (|ST|) in 

S = {{x, t) : -r < t < 0, /i + - 5 < \x\ < h + ut + r/}, 

where r] > and 5 > small are to be chosen and 

h — h T + r] 



T 



so that will cross it at Pi, yielding a contradiction. 



Consider a spatially radially symmetric function 



(p(x,t) = f(\x\ - ujt) 



f(s) 



2c 



(-Ks + a) 3/2 + b. 



3K 



where 




and the constants a, b £ R and c > are to be chosen below. 



Observe that g(s) = (f'(s)/c) 2 = —Ks + a. Choose < i] < |/i r small enough and a € R 
so that \ < g(s) < 1 for all h — r] < s < h + rj. It is also clear that we can find c > and 
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b £ M such that < f[h) < L and f(h + rj) > M. / is continuous and so there is also 5, 
< 5 < r) such that f(h - 5) > 0. 

In particular, with this choice of i], 5, a, b and c, /(s) is positive and strictly increasing in 
[h — 8, h + ?/] . It follows that (/> satisfies 



To conclude, we note that 4> > u on the parabolic boundary of S. This is clear from the fact 
that (f> > in S and cf) > M on those parts of the parabolic boundary of S that intersect 
{u > 0}. But (u — (f))(Pi) > and therefore u crosses (j) from below inside S. That is a 
contradiction with u being a viscosity subsolution of (|ST|) . 

If ti > to) including the case when (P\,Pq) is vertical, the construction of a supersolution 
is more straightforward but we can also use the above construction with h T = h. □ 

Corollary 3.7. The vectors Pq — Pi and Po — P2 cannot be vertical (i.e., be parallel to the 
time- direction). In other words, the space projections of P± and P2 cannot be equal to the 
space projection of Pq . 

Proof. Observe that if say t\ < to and the line (Pi, Po) is vertical, then also (P2, Po) must be 
vertical and t<i > to. This follows from the strict ordering Z(-,t) < W(-,t) for all t < to and 
B\ being an interior ball of {Z > 0} and B2 being an interior ball of {W < 0}, i.e. Z > 
in Pi and W < in B 2 which implies Pi n P 2 n {t < t } = 0. That is, Pi n P 2 = {Po} 
whenever one of the balls is contained in {t < to}. But Lemma 13.51 asserts that (P2,Po) 
cannot be vertical with P 2 in the future of Pq. Similarly when (P 2 ,Po) is vertical with 



Due to Corollary 13 .7\ the interior ball Pi of {Z > 0} has an interior normal (ei,m) at 
Po £ dB\ with \m\ < 00. Note that m denotes the (outward) normal velocity of the space- 
time ball Pi at Po. Therefore, formally speaking, m is smaller than the normal velocity V n 
of {Z > 0} at Po. Since Z is a subsolution of (|ST|) and since Po £ d{Z > 0} = d{Z < 0} , 
from the discussion in section 2 (see (V)) formally it must be true that 




> 



since to + ^-4- < -y in S. Therefore <j> is positive and strictly supercaloric in S. 



t 2 < t . 



□ 



m < \DZ + \ - \DZ~\. 
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Next lemma states that this is indeed true in a weak sense. 



Lemma 3.8. 

. t Z{P + \ ei ,0) . Z(P -/xei,0) ^ 

lirnmf h limmt > m. 

A->0+ A fi^O+ fl 

Proof. Let us define 

. Z(P + Aei,0) . Z(P -/xei,0) 

a := limmt ana o:=nmmt . 

A->0+ A n->o+ n 

Observe that b < 0. Indeed by Corollary 13 .7\ we can find a smaller closed ball B G P>2 such 
that 5 n <9P 2 = Po and /i > such that Bn{t = t -li} / 0. Since Z < W for t < t , 
there is a smooth function g > that solves the heat equation inside P n {to — h <t < to} 
such that Z + 5 < W on B f) {t = t — h}, g > in B° C\ {t = t - h} and g = on 
<9P n {to — ^ < t < to}. As Z is subcaloric and is supercaloric in B n {to — /1 < t < to} 
it follows that Z + g<WonPn{t = to} and 

> hminf ^ (P °-^ ei ' 0) > b + \Dg (P )| > 6. 

Suppose a + b < m. Since & < 0, e > can be chosen sufficiently small so that 

a + 6 < m — 6e and 6 + 3e < 0. 

We will construct a supersolution of (|ST|) constructed in a parabolic neighborhood K of 
Pi, which will be specified later, which crosses U from above, to yield a contradiction. We 
point out that, due to the two-phase nature of the problem, the construction of such barrier 
and the domain K requires modifications from the one-phase version constructed in [10] 
and [I?]. 

From the definitions of a and b, there are sequences A& — > 0, Hk — > with A& > and > 
such that 

Z(x + A fc ei, t ) < (a + e)A fc , 
Z(x - jt/fcei, t ) < (6 + e)/u fc . 

In what follows, A £ {A/ c }^ =1 and /i G {/ifc}j£Li- Now denote Q = Po + (Aei,0), P = 
Po + (— /iei,0) and Pq, resp. Pr the balls centered at Q, resp. P, with radius r. The 
situation is illustrated in Figure El 

From the definition of Z : it follows that U(x,t) < Z(y,s) for any (y,s) G B r (x,t) and thus 

'Z(Q)<(a + e)A in P Q , 
P < < Z(P ) = in P , 

Z{R)<(b + e)n in B R . 
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Figure 5: The slice {t = t±} 



In fact, U < on a slightly larger set. Because Z < W at t = to and W < V(I%) = in B2, 
we can conclude that Z < on D = B2 n {t = to}- It means that £/ < on E, the union of 
balls of radius r with a center in D, 



£ = |J B P (z,i) 

(x,t)S-D 



Finally, we have an improved bound, 



f/ < < 



(a + e)A in Sq, 
in E, 

K (b + e)fi in B R . 



(3.5) 



Now we will compare U with a supersolution in a domain := Kr . Fix r > small and 
define 

Ky = {h-r<t<t 1 }n \J (SB \iiitE-7(ei,0))U5o\intJ3fl. 

0<7<A 

In dimension n = 2, the shape of this set is depicted in Figure [6] and the shape of its slice 
{t = ti} is depicted as the thick curve in Figure[5j The boundary of the slice is given by the 
arc S = 8Bq \ E, the arc T = 8Br n Bq, and finally of the curve (surface) C that connects 
S and T. C at any given time is composed of C\ = dL n E \ Br, where L is a cylinder with 
axis ei connecting / = 8Bq n <9E with (9i?0) and C2 = 8Bq \ (Br U L). 

Let <f>(x,t) = <fi(p,t), where p = \x — xq\. We need U < (j) on the parabolic boundary of 
Kr ,fl . By Taylor expansion in A and (i — it follows that the spatial distance of the point 
of intersection / from xq can be found explicitly for each time t as 

Rj(t) = di + A - m(t -h) + o (A) + o(t - h) < di + X - m(t - h), 
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Figure 6: Domain K for n = 2 

with d\ from (|3.4p . Moreover, Rxit) < \ x ~ x o\ < Rs(t) m where 

= di + A - m(i - h) + o(i - h) and i? T (0 = di - /x - m(i - ti) + o(i - t x ). 

Lastly, the radius of the space-ball Bq n {t = t} is given by 

R^t) =di- m(t - h) + o(t - ti). 

Now we are ready for the construction of the supersolution. Take a = a + 3e, (3 = b + 3e, 
R = di = Ri(t\) and m = a + /3 < m. Let 0(x, t) be the function <fi°~ _ (a; — xo, i — io) from 
Appendix. Let t > 0, 5 > be the small constants and S T ( 5 be the set from Appendix 
where <fi exists and is a radially symmetric classical supersolution that is "e-close" (in the 
sense of (14.5|) ) to linear functions with slopes a and (3 on each side of the spherical free 
boundary of <fi moving with speed fh and with radius R at t = t\. 

Let R(t) = R — rhit — t\) be the radius of the free boundary of (j>. We claim that for 
sufficiently small A, [i and r we have U < 4> on the parabolic boundary of K T '^. 

To prove the claim, first note that for sufficiently small A, fj, > 

c S Tj5 = {(s,t) : R(t) -5 <\x-x \< R(t) + 5, h - r < t < ti}. 

We will compare <j) with U on two parts of the parabolic boundary. 
Part 1: On lateral boundary d L (Kr'^ D {h - r < t < ti}): 
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Prom (|3.5p we need 



(a + e)A for R r (t) < p < Rs(t), 
0(p,t)>lo for Ri(t) < p < Rj(t), 

\b + e)p for R T (t) < p < Pi(i). 



(3.6) 



The second condition is satisfied as long as R(t) < R\(t) and that is true for all t £ [t\—T, t\] 
if r is small enough. Third condition is trivial for R(t) < p < Ri(t). For Rx{t) < p < R(t) 
we use the fact that 



<f>{p,t) >(/?-£) [R(t)-p 

>(/3-e) (R(t)-R T (t) 
> (P-e)p > (b + e)p. 

Lastly the first condition in (|3.6p holds since 



Hp, t)>(a- e) (p - R(t)j >{a-e) (P/(i) - R{t) ^ 

> (a — e)(d\ — m(t — t\) + A — r\ — d\ + rh(t — ti)) 
= (a — e) ((m — m)(t — t\) + A — rj) 

> (a - e)(A - ??) > (a + e)A 

if r and A are small compared to e since 77 = o(A) + o(t — t\). 

Part 2: On r\{t = t 1 -r}: 

From the bounds on U in (13.51) we need 



<P(p,t) > 



(a + e)\ for R^tx - r) < p < R s (ti - r), 
for P T (*i -r) < p< -r). 



The second inequality can be satisfied by taking r and p small so that R(ti — r) < Rt^i — t). 
The first inequality is then satisfied by taking A small. 

We demonstrated that U < 4> on the parabolic boundary of Kr ,fl and U > at Pi if A, 
and t > are small enough. Since is in fact a strict classical supersolution due to 
Appendix, we can perturb (ft into 4> near Pi such that (p is still a classical supersolution with 
U < <ft on the parabolic boundary of i^r'^ and U > at Pi. This is a contradiction with 
the definition of U. 

□ 
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We recall that Z < W at t = to- Therefore Lemma 13.81 as well as considering g in the 
beginning of its proof yields 

W(P + \ ei ,0) W(P -» ei ,0) 
lim sup h nm sup (3. 7 J 

> Ummf Z(Po + A ei ,0) + 1 , m , nfZ (Po-Me 1 ,0) + 

> TO. 



Recall that strict ordering Z <W for t < to implies that B<i has an exterior normal (ei, m) 
with to < to. With this result in mind, we can repeat the argument in Lemma 13.81 to get a 
contradiction by proving that the quantity (|3.7p must be less than or equal to m. 

This establishes the desired contradiction and we conclude that Z and W stay strictly 
separated for all times. 

Note that Theorem 13. II yields ordering only between strictly separated initial data. There- 
fore it does not directly yield the uniqueness result for viscosity solutions, even with the 
additional assumption \{uo = 0}| =0. However under some restrictions on initial data, 
uniqueness readily follows from Theorem 13.11 

Corollary 3.9. Letu(x,t) be a viscosity solution of ()STp with initial data uq and boundary 
data 9(x,t) = 9{x) < 0, where uq and 9 satisfies one of the following conditions: 

(a) uq is strictly concave; 

(b) 9(x,t) = 9 < and both f2 and uo is star-shaped: i.e. for sufficiently small e > 0, 

(1 + e) _1 C VL and uq(x) < uq((1 + e)~ 1 x), where equality only holds at x = 0. 



Then u(x, t) is unique. 



Proof. 1. Let us begin by mentioning that for any a, b and c > 0, if u is a viscosity sub 
(super)-solution of ()ST|) . then so is 

au(bx, ab 2 t + c). 



2. Now let u and v be two viscosity solutions of ()STp with initial data uq. In the case of 
(a), one can verify via barrier arguments that v(x,t) < uq(x) if < s < t with equality 
holding only at 8lQ- In particular, for sufficiently small e and 5 such that 0<e<<5cl, 
we have (1 + e)v(x,5) < uq(x) (recall that 9 is negative). Therefore observation made in 
step 1 and Theorem 13.11 yields 

(1 + e)v(x, (1 + s)t + d) < u(x,t) in Q 
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for any < e <C <5 <C 1. Now we can send 6 — > and e — > to obtain v < u. Similarly we 
can argue with t> and -u switched to obtain u < v, and therefore u = v. 

In the second case, note that for any e > and 5 > 0, 

(1 + (fT^oOc) < u ((l + in CI. 

Also, since tio is star-shaped, the minimum of uq (and u and f ) equals 9. Since u and v are 
negative caloric functions near dfl, Hopf's lemma yields that the inward normal derivative 
of u(-,t) and v(-,t) on d£l is strictly positive. In particular for any T > there exists 
5 = 0{e) such that 

(1 + 5)~H < u((l + e)~ l x, t) for (x, i) e d L Q. 

Hence Theorem 13.11 again yields that, for sufficiently small e > 0, 

(1 + 8)- x v{x, (1 + < + e)- 1 ^, (1 + e)" 2 t) in Q. 

Therefore, by sending e — )• 0, we can conclude that v < u. As before, by switching v and u 
in above argument one can conclude v = u. 

□ 

For a general uniqueness result using the comparison principle we refer to section H] in |17j 
where rather technical barrier arguments take place. We have not attempted to follow 
this approach since, as mentioned before, we are interested in the coincidence of the two 
generalized solutions for ()STp . 

4 Correspondence between weak and viscosity solutions 

Recall that, due to [I5j[18j[19] > there exists a unique weak solution h of problem (jSTp . defined 
in Definition 12. 1^ with initial data ho and boundary data 6. Let us define u = x(h) to be 
the weak temperature solution of (jSTp associated with initial enthalpy function ho, see |15j . 

Lemma 4.1. Any weak temperature solution u of (jSTp with initial data uo = x(ho) is a 
viscosity solution of (|ST|) . 

Proof. The analysis in [8] yields that weak temperature solutions are continuous for contin- 
uous x(^-o)- Moreover, a comparison principle holds for weak solutions. Indeed, if h\ and 
/i2 are two weak solutions with h\ < hi on dpQ, then h\ < hi in Q. See [19], for instance. 

The comparison principle can be extended to accommodate weak sub and supersolutions. 
A function hi is a weak subsolution if it satisfies (|2.ip in Definition 12.11 with > instead of 
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2 1 

equality for all tp S W 2 ' (Q), p > 0, p = on 8lQ U {t = T}. A function h,2 is a weak 
super solution if it satisfies (|2.ip with < instead of equality for all <p as above. A comparison 
principle still holds. If h% < hi on the parabolic boundary dpQ then hi < h 2 in Q. It is 
straightforward to check that the comparison principle on a cylindrical domain implies a 
comparison principle on a general parabolic neighborhood in the sense of Definition 12.21 

A standard integration by parts also shows that a classical subsolution 4> in Definition 12.31 
corresponds to a weak subsolution x~ 1 ( ( / , )i an d also that a classical supersolution tp is a 
weak supersolution note that, by definition, <f> and p have no mushy region and 

thus x~ 1( P an( i X f are wen defined almost everywhere. We deduce from this that we 
can compare weak temperature solutions with classical sub and supersolutions. Hence we 
conclude that weak temperature solutions are viscosity solutions of (|STj) in the sense of 
Definition [231 □ 

Corollary 4.2. There exists a viscosity solution of (|ST|) for continuous x(ho) and 9. 



Unfortunately it is not clear to the authors how to find a corresponding weak solution for 
a given viscosity solution when the initial data has a fat zero set. However it is possible 
to identify the maximal and minimal viscosity solution with maximal and minimal weak 
temperature solutions (see Proposition 14. 4p . To prove this we first need a uniform control 
on the expansion of the sets {u < 0} and {u > 0} at t = in terms of the L°°-norm of the 
initial data. 



Lemma 4.3. Let u be a viscosity solution of ()STp with initial data uq. Then for any 
r > the set {u < 0} does not expand by more than r before the time t r = cr 2 , where 
c = c(n, M) is a constant specified later in the proof, depending only on the dimension n 
and M = maxg p Q \uq\. In other words, 

{x : u(x, t) < 0} C {x : dist(x, {u < 0}) < r} 

for all t G [0, t r ] . 

The same is true for {u > 0} . 

Proof We only prove the lemma for {u < 0}. Parallel argument holds for {u > 0}. 
Set 

A r = {x : dist(a;, {uo < 0}) > r}. 
We need to show that u(x,t) > in A r x [0, t r ], for t r specified below. 
Thus choose xq £ A r . Define the function 

jp-t - 1, n > 3, 



p" 

— log p, n = 2. 
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Note that / is decreasing and /(l) = 0. For simplicity we show the following computation 
for n > 3 only. The case n = 2 is similar. 

Set K = hMji and v = 8K (™~ 2 > _ "w e w [\\ construct a classical subsolution of (|ST|) in the set 

C = {(x, t) : \x — xq\ < r — vt, < t < t r } 
that is below n on the parabolic boundary of C and is positive at xo for < t < t r . 
Define the radially symmetric function <f> = (f)(x,t) in C as 



<f>(x,t) 



'Kf fel), 2|x-x |>r-^ 



r—vt 

x, t), 2\x — xq\ < r — vt, 



where 4>(x, t) is a solution of the heat equation in the domain 

Ci = | i x , t) : \x - xq\ < i (r - vt) , i G (0, t r ) 

with smooth positive initial data 4>(x, 0) < u(x, 0) in the interior of D = {x : \x — xq\ < r /2} 
and boundary data 0. 

The modulus of the gradient |-D</> _ | on the free boundary T((f>) is 

\D*-\-K\rm\ 2 - 2K(n - 2) 

For t < t r = £ = c(n, M)r 2 , 



vt r — vt 



, n ,_, 4*T(n-2) u 
ii r 2 



Therefore is a classical subsolution of (IST|) in the set C. Also, < n on the parabolic 
boundary dpC. The definition of viscosity subsolution then implies <f) < u in the set C and 
therefore it > at xo for < t < t r . □ 

Proposition 4.4. There exist maximal and minimal viscosity solutions U\(x, t) and U2(x, t) 
of (ISTj) satisfying the following: 



U2 < u < C7i /or any viscosity solution u of (jSTp . 
Furthermore, U\ is the unique weak temperature solution of (jSTj) luii/t i/ie initial enthalpy 



^ max . _ 

n 



uo if u > 0, 

n - 1 if n < 0, 



and t/2 ^ e unique weak temperature solution of (jSTp u>ii/i i/ie initial enthalpy 



trmra 

n 



n £/ n > 0, 

n - 1 if n < 0. 
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Proof. 1. Take an approximation Wq < uq < Vq such that the zero set of the approximation 
has zero measure. This can be done, for example, with simply Wq = uq — e and Vq = vq + e 
for along a sequence of e — > 0. Indeed the function F(e) := \{u > 0}| — |{uo > e}\ is 
a monotone, bounded function, therefore it can only have countable number of jumps. 
Except at these jumps, |{uo = e}\ = 0. 

Let us choose a viscosity solution v £ (and w £ ) of ()STp respectively with initial data Vq 
(resp. Wq) and lateral boundary data 8 + e (resp. 6 — e). In particular (excluding a 
countable number of e) the functions Vq and have no mushy region, and therefore h^' 6 := 
X _1 (^o) ana - ^o' £ := X~ 1 ( v o) are well-defined almost everywhere, hence the corresponding 
weak solutions h 1,e and h?' £ exist. Now, by Lemma 14.11 the weak temperature solution 
v £ := xih 1 ' 6 ) (resp. w e := x(^ 2 ' £ )) is a viscosity solution of (jSTp with initial data Vq (resp. 
Wq). Let u be a viscosity solution of ()STp with initial data uo- Let w be the unique weak 
temperature solution of (IST|) with boundary data /io := X~ l { u o)- For convenience, let us 
denote the initial and lateral boundary data g := (ito, ^). 

2. The proof of the following lemma will be presented in the Appendix: 

Lemma 4.5. The functions v £ (orw e ), along a subsequence, uniformly converge to a con- 
tinuous function u with boundary data g . 

Due to the stability of viscosity solutions, it is straightforward to verify that u is a viscosity 
solution of (1ST!) with boundary data g. 

On the other hand, due to stability result of [19], for each time t > the weak soution 
h r ,£ (-,t) converges to h(-,t) in L 1 -norm, where h(-,0) = h™^ . 

Using the uniform convergence of v £ := x(h 1,£ ), it is easy to check that u = x(h) = U\ a.e. 

3. Lastly, Theorem 13.11 yields that any viscosity solution u of (jSTft with initial data uq 
satisfies 

w £ < u < v £ , and thus taking e — > we obtain U2 < u < U\. 

□ 

The following Theorem is obtained as a corollary of Propostion 14.41 

Theorem 4.6. The viscosity solution u of ()STp with \{u$ = 0}| =0 coincides with the 
weak temperature solution and it is unique. 
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Figure 7: Test function <f){r, t) at time t 

Appendix 



First we construct "two-phase" test functions for the barrier arguments in the proof of the 
comparison principle, Theorem 13.11 

A . Test functions. Here we provide a construction of a class of radially symmetric classical 
sub- and supersolutions for that purpose. We will be working in spherical coordinates and 
thus we shall denote p = \x\. 

Let a / 0, e > 0, R > and m be given parameters satisfying 

sign a ^ sign/3 and a + f3 ^ ?nsign/3. (4-1) 

We show below that one can find constants r > 0, 5 > and q such that a test function 
(j>(x, t) of the form 

4>(x, t) = 4>o{\x\ — R — mt). 



with 



(s) 



as + q^- s > 0, 
-f3s + q^ s < 0, 



defined in the domain 

S ri5 = {{x,t) : R + mt - 5 < \x\ < R + mt + 5, -r < t < 0}, 

is a classical subsolution (or supersolution) of (|ST|) on S Tj( 5, see Figure [71 Moreover, 4>(x,t) 
is "e-close" to a linear function of p on both sides of the free boundary T(0), see (|4.5p . 

We shall denote Rt = R + mt since the set T(cp) = {(x, t) : r = R + mt} is the free boundary 
of (p. Observe that on T(0), we have 

I D(j) + 1 = max (a, /3) , I D(fi~ I = — min (a, /3) . 
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and thus 



\Dcp + \ - \D<j>~\ =a + p. 



(4.2) 



We can also deduce the direction of the outer normal of {(f) > 0}, and hence the sign of 
the normal velocity V n , from the signs of a and /3 due to (|4. 1 1) . As T((p) is the space ball 
expanding with speed m, we recognize two situations: 



Vn 



—m when (3 < 0, 
m when /3 > 0. 



(4.3) 



Now we can express the free boundary condition in the light of (|4.2|) and (|4.3|) as 

(a) m sign (3 > a + (3 for a classical supersolution, 

(b) m sign /3 < a + (3 for a classical subsolution. 

Suppose that the condition in (a) is satisfied and we are thus constructing a classical su- 
persolution. In this case, the function <j> has to satisfy also (dt — A)</> > in {<f> ^ 0}. 



n 



1 



Rt 

n — 1 



Straightforward computation yields 

{d t - A)(f>(x,t)\ p=R + = -a | 

(d t -A)cf>(x,t)\ p=R -=[3\m + 

Continuity allows us to choose small r > and 5 > and a suitable constant q such that 

in E T ,a \ r(^) (4.4) 



and 



(d t - A) 4>{x, t) > 

/8 



4>(x,t) 



< e. 



Rt 



< e, 



when p > R t , 
when p < R t , 



(4.5) 



in S r a. 



Thus let 4>^ m denote the classical supersolution constructed above with the parameters R, 
m, a and [3. Also let 5 [4 >C Rm> e l an< ^ T (^??'m' e ) denote the constants from above. 



The case when the condition in (b) is satisfied and we are constructing a subsolution is 
analogous. The constants r > 0, 5 > and q are chosen in such a way that the function 
4> satisfies (dt — A)cj)(x,t) < in (I4.4p . And correspondingly, let 4>R m ' 3e * ne constructed 
classical subsolution. □ 
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We next present the proof of the lemma used in the proof of Proposition 14.41 

B. Proof of Lemma \4-5\ We will only prove the Lemma for v £ . The argument for w £ is 
parallel. Let u k := v l ' k and let g k the corresponding initial and lateral boundary data of 

[8] provide us with a uniform bound on the W^iQ) norm and with a uniform modulus of 
continuity for v e in the domains away from dpQ: 

Q d = {(x, t)€Q: dist((x, t),d P Q) >d}, d> 0. 

Therefore we need a uniform modulus of continuity on dpQ to be able to extract a subse- 
quence of u k converging uniformly. We will obtain this by comparing u k to solutions of the 
heat equation near the parabolic boundary of Q. 

First find ko for which g k ° < on dj_,Q. Theorem 13.11 implies u k decreases in k. In particular 

|n fc+1 > O} C \u k > O} C \u k ° > o}. 
Since {u k ° > 0} is a closed set and g k < for k > ko, there exists rj > such that 
> J C Q TI := {(x, t) G Q : dist(x, dtt) > r]} for all k > k . 

Note that Q" 1 := x (0,T) where 

fi' = {i£fl: dist(x,aO) > ??}. 

Define ^i, resp. to be the solution of the heat equation in Q with boundary data 
resp. —g k ~- Here <7 fc+ and g k ~ are the positive part and the negative part of g k , 
respectively. 

The maximum principle for the heat equation yields 

■\\) k _ < u k < in Q. (4.6) 

A uniform modulus of continuity for u k on the lateral boundary OlQ then follows from 
comparing u k with tp k that solves the heat equation 

'<pt-A<p = 0, inQ\Q J? , 
< cp = g k , on d L Q, 

K <p(x,0) = g k (x,0)((x), onQ\Q T '. 

Here C( x ) is a smooth function on $7 \ f]* 7 such that C G [0, 1], ^ = 1 on \ ri 77 / 2 and £ = 
on dW. 
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Observe that (|4.6p and comparison for the heat equation implies that 

^-<u k <(f k on Q\Q T >. 
Since ip^_ = u k = ip k = g k on 8lQ, the uniform continuity follows. 

The proof for t = is more involved. We split u k = u k+ — u k ~ and consider each of u k± > 
separately. 

Choose r > small and fix k. Let A k = {x G $7 : g k > 0}. Define 

A k = {x e A: dist(x, dA n Q) > r}. 

By Lemma |4.3| there is i r > such that A k C : u fc (x, £) > 0} for all t < i r . Moreover, i r 
can be chosen independent of k. Construct a solution ip k to the heat equation in A k x (0, t r ) 
with boundary data on (dA k n J7) x (0,t r ) and g k (x,t)C k {x) on the rest of the parabolic 
boundary of A k x (0, t r ). C, k is chosen to be a smooth cutoff function with values in [0, 1] 
on A k such that = 1 on A k r and ( fc = on dA k flJJ. It is also chosen in such a way that 
||C fc ||co,i is uniformly bounded in k. Define ip k = in the complement of A k x (0, t r ). The 
maximum principle for heat equation then yields 

<Pr < u k+ < ^+ in O x [0,< r ]- 
Thus the standard parabolic estimates for the heat equation (|18j) imply 

u k+ { Xl ,h) -n fc+ (x ,0)| = |^ + (xi,ti)- 5 fc+ (x ) 

< w 9 (2r) + (jj H (\xi - x \ 2 + hj , 

for all t\ < t r , xo,x± S fi, modulus of continuity wjj independent of k and r, where u g is a 
modulus of continuity of g. This implies uniform modulus of continuity for u k+ at t = 0. 

A similar consideration yields a uniform modulus for at t = 0. Those together give us 
a uniform modulus of continuity for u k . 

Now as in [8], we can find a subsequence {u kj }°° =l such that 

f3 (V^ -> /3 («) weakly in L 2 (Q), 

Du fcj ' -)■ Du weakly in L 2 (Q), 

u kj — > u uniformly on Q. 

□ 
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